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We present a method of construction of exact localized many-body eigenstates of the Hubbard 
model in decorated lattices, both for U = 0 and U —>■ oo. These states are localized in what concerns 
both hole and particle movement. The starting point of the method is the construction of a plaquette 
or a set of plaquettes with a higher symmetry than that of the whole lattice. Using a simple set 
of rules, the tight-binding localized state in such a plaquette can be divided, folded and unfolded 
to new plaquette geometries. This set of rules is also valid for the construction of a localized state 
for one hole in the 17 —^ oo limit of the same plaquette, assuming a spin configuration which is a 
uniform linear combination of all possible permutations of the set of spins in the plaquette. 
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The field of itinerant geometrically frustrated elec¬ 
tronic systems has attracted considerable interest in the 
last two decade^i^f^. Much of this interest was re¬ 
lated with the study of fiat-band ferromagnetism in these 
system^i^ll^. Flat-band ferromagnetism occurs in deco¬ 
rated lattices of the Mielke’s and Tasaki’s classes, which 
display degenerate localized CTOund states with overlap¬ 
ping probability densitie^^mSS) emerging ferromag¬ 
netism can be interpreted as resulting from a general¬ 
ized Hund’s rul^^. In the case of the lattices which fall 
into the Lieb’s class, the fiat bands intercalate itinerant 
bandd^ and mean-field studies of the Hubbard Hamil¬ 
tonian in the Lieb lattice indicate that for large U, fer¬ 
romagnetism is expected ex cept near half-filling where a 
ferrimagnetic phase appears^SESl 

These localized states are one-particle eigenstates of 
the tight-binding Hamiltonians for the decorated lattices 
and little is known about the many-body eigenstates of an 
interacting system of fermions in decorated lattice^^ (as¬ 
suming Hubbard-like interactions), besides the appear¬ 
ance of a ferromagnetic ground state i n de corated lat¬ 
tices of the Mielke’s and Tasaki’s classe^^*22) particu¬ 
lar, the interacting ground state of the Hubbard model is 
not known in the case of lattices of the Lieb’s class. Ap¬ 
proximate analytic results can be obtained in the weak 
coupling limit, addressing interactions as perturbations 
of the tight-binding fiat bandd^. 

In this manuscript, we present a method of construc¬ 
tion of exact localized many-body eigenstates of the Hub¬ 
bard model in decorated lattices of arbitrary dimensions, 
for [7 = 0 and U —>■ oo. These states are localized in 
what concerns hole and particle movement. This method 
relies in simple arguments which lead to a set of quan¬ 
tum “origami” rules: i) if one plaquette or a set of pla¬ 
quettes has a higher symmetry than that of the whole 
lattice, one can find energy eigenstates that have zero 
probability density at the sites that connect the plaque¬ 
tte or the set of plaquettes to the rest of the lattice (this 
argument is enough to justify the existence of localized 


states in the case of two-dimensional decorated lattices 
of the Lieb’s class); ii) given such a localized state in the 
symmetric plaquette, one can fold the plaquette, either 
at the probability density nodes or at other equivalent 
sites (adjusting the probability density at those sites and 
the hopping constants that involve those sites), therefore 
lowering the symmetry of the plaquette; iii) the energy 
of the localized state can be lowered by adding hopping 
terms between sites with the same localized state phase 
(if the hopping constant is negative) or hopping terms 
between sites with opposite phases (if the hopping con¬ 
stant is positive). Hopping terms between nodes of the 
localized state may also be added, but do not change the 
energy of the localized state. The hopping terms added 
must preserve the symmetry of the localized state. These 
two arguments justify localized states in decorated lat¬ 
tices of the Mielke’s and Tasaki’s classes. Furthermore, 
the spin degree of freedom of the U = 0 Hubbard Hamil¬ 
tonian may be interpreted as a sublattice index and local¬ 
ized states can also be created using these origami rules 
involving the two (up and down spin) sublattices. Such 
localized states arise for instance as edge states in ID 
tight-binding descriptions of topological insulator^^MH! 

The remaining part of this paper is organized in the 
following way. First, we review the construction of one- 
particle localized eigenstates of the tight-binding deco¬ 
rated lattices of the Lieb’s class. We then generalize this 
construction to more complex lattices using a symme¬ 
try argument and introducing the set of origami rules. 
Next, we show how to extend these rules to the case of 
the U —i' oo limit of the Hubbard model. Finally we 
conclude. 

The Hubbard Hamiltonian in a decorated lattice can 
be written as 

H= 'y ^ + 17 (1) 

(bRo- i 

where the creation (annihilation) of an electron at site i 
with spin a is denoted by {cia) with riia being the 
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number operator riia = c|^Cicr and = rii^ + The 
sum over {ij) is the sum over all pairs of sites with a finite 
hopping probability between them and this a different 
sum for each decorated lattice. The hopping constants 
are assumed to be equal, = t, unless stated otherwise. 
When t = 0, all states with the same number Nd of 
doubly occupied sites are degenerate. In this paper, we 
assume Nd = 0. The Hubbard model in the limit U —>■ oo 
is also designated as Harris-Lange modeP^. In this limit, 
using the identity Cia- = Cicr[(l — nia-) + riia-], the Hubbard 
model can be rewritten as 

H 'y ^ T t/(2) 

(b'l.o’ * 

with a = —a. An important point about the strong 
coupling limit is that the Hamiltonian eigenfunctions, in 
the case of a Hubbard ring, can be written as a tensorial 
product of the eigenfunctions of a tight-binding model of 
independent spinless fermions (holes) in the ring with L 
sites and the eigenfunctions of an Heisenberg model (with 
exchange constant J = /U) in a reduced chaicpSHSI. 

Let us first discuss the U = 0 case. Flat bands in the 
one-particle tight-binding energy dispersion of geomet¬ 
rically frustrated lattices reflect the existence of degen¬ 
erate localized eigenstates which are translated versions 
of the same state IV^ioc)- The probability density associ¬ 
ated with one of these localized states is non-zero only 
in a small lattice region. In the particular case of deco¬ 
rated lattices of the Lieb’s class, the localized states can 
be viewed as one-dimensional standing waves in tight- 
binding rings, associated with paths in the 2D lattice 
which include one or two plaquette^^. For zero flux, all 
one-particle energy levels of a tight-binding ring (with 
even number of sites) are doubly degenerate (except for 
k = 0 and k = ir) and the respective eigenstates have 
opposite momenta. Adding or subtracting the states of 
opposite momenta, one obtains a standing wave with a 
number of nodes that depends on k. If these nodes co¬ 
incide with the sites at the vertices of a plaquette of a 
decorated lattice, the electron becomes trapped inside 
the plaquette. Therefore, flat band eigenstates of dec¬ 
orated lattices of the Lieb’s class are constructed from 
standing waves such that the nodes coincide with sites at 
the vertices. Note that these localized states overlap in 
real space, that is, they constitute a basis of the subspace 
of localized states but not an orthogonal basis. 

The previous argument for lattices for the Lieb’s class 
can be generalized to decorated lattices of the Mielke’s 
and Tasaki’s classes and other decorated lattices using a 
symmetry argument. First, let us discuss the case of the 
Lieb lattice (see Fig. [^. The tight-binding Hamiltonian 
of one plaquette of the Lieb lattice has the symmetry of 
a ring of 8 sites, that is, the plaquette Hamiltonian is 
invariant in a 2tt/S rotation of the set of site indices (or 
equivalently in a circular permutation of the set of site 
indices). We emphasize that this rotation should not be 
confused with a 2tt/8 rotation in real space (the plaquette 
is not invariant in such a rotation). However, the rotation 




Figure 1. The symmetry of the tight-binding Hamiltonian for 
the Lieb plaquette is the same as that of the tight-binding 
ring and larger than that of the Lieb lattice. The distance 
between adjacent sites is assumed to be a = 1. All hopping 
constants are equal. 


of 27r/4 in the set of sites indices can be interpreted as a 
27r/4 rotation in real space. In the ring of 8 sites, one has 
outer sites (that are connected to the rest of the lattice) 
and inner sites (with connections only to sites of the pla¬ 
quette). In the case of Fig. the inner sites are indicated 
by the filled circles and the outer sites are given by the 
empty circles. The generator of this rotation symmetry is 
the equivalent of the angular momentum in the ring (note 
that in a 2D lattice, the direction of the angular momen¬ 
tum is always perpendicular to the lattice and therefore 
equal to mh, where m can be interpreted as the m in the 
ring momentum k = m ■ 2 'k/N) and one can construct 
an eigenbasis of the Hamiltonian which is simultaneously 
an eigenbasis of the angular momentum. The time re¬ 
versal symmetry of the Hamiltonian implies that each 
eigenstate of the plaquette tight-binding Hamiltonian is 
degenerate with the respective state obtained in a time 
reversal and these states have opposite angular momenta 
(this is equivalent to stating that ring eigenstates with 
momenta k and —k are degenerate). These two states 
can be added or subtracted, generating the equivalent 
of the standing waves in the ring, that is, states with 
zero probability density at certain sites of the cluster. If 
the angular momentum is HN/A, where N is the number 
of sites of the ring, one has zero probability density at 
the inner sites or at the outer sites of the Lieb plaque¬ 
tte (nodes are separated by A/2 = 2). The latter will 
be a localized eigenstate not only of the Lieb plaquette 
but also of the tight-binding Hamiltonian of the full lat¬ 
tice. Note that this description is valid for any plaquettes 
which have the same rotation symmetry as the ring. For 
example, one could add additional sites at the center of 
the Lieb plaquette and the rotation symmetry would re¬ 
main, as shown in Fig. (in all Figures, the relative size 
of the circles that represent lattice sites corresponds to 
the relative value of the wavefunction amplitudes on the 
sites). 

Thus, our first rule is that localized states can be con¬ 
structed if a plaquette (or a set of adjacent plaquettes) 
has a larger symmetry than the lattice, so that the Hamil¬ 
tonian has two degenerate eigenstates (which are simul- 
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Figure 2. In these simple variations that retain the symmetry 
of the Lieb plaquette, two localized states exist corresponding 
to two possible choices of angular momentum. The size of the 
circles that represent lattice sites indicates the relative value 
of the wavefunction amplitudes. There is an analogy of these 
states with atomic orbitals. 


taneously eigenstates of the generators of the symmetry 
of the lattice) which have the same wavefunction val¬ 
ues at the outer sites of the plaquette (see Fig. Sal. 
This rule is enough to explain the existence of localized 
states in lattices of the Lieb’s class. More complex lat¬ 
tices with localized states can be constructed by adding 
sites or hopping bonds that do not lower the symmetry 
of the plaquette. These additional hoppings can be di¬ 
vided into two sets: i) the set of hoppings from or to 
sites with probability density nodes (these hoppings do 
not modify the energy of the localized state); ii) the set 
of hoppings between sites with finite density probability 
(these hoppings lower or raise the energy of the localized 
state). 

A second rule for the construction of lattices with lo¬ 
calized states is the following. The existence of sites 
where a localized state has probability density nodes does 
not affect the energy of the state and these sites can be 
dropped, duplicated (as well as the respective hopping 
bonds), or simply added (introducing appropriate hop¬ 
pings with neighboring sites) and the localized state re¬ 
mains an eigenstate of the modified tight-binding model 
associated with the new plaquette geometry (see Fig. 


3b I. Furthermore, if one can draw an axis through the 


plaquette that crosses only nodes, then dropping these 
nodes one divides the localized state into two eigenstates 
of the tight-binding Hamiltonians associated with the 
parts of the plaquette. Bonds between nodes can also 
be dropped, added or duplicated. This rule justifies the 
localized states in the lattice of Fig. 4a In fact, sites 
A and B in Fig. can be seen as a duplication of the 
equivalent site of the ring of Fig. with the addition of 
a hopping bond between the duplicated nodes. 

The third rule consists of the following: localized states 
can be folded (adjusting the amplitude at the crossing 
and the respective hoppings) along an axis that crosses 
the plaquette through sites that have the same wavefunc¬ 
tion values (see Fig. 3cI; if the folding is along an axis 
that crosses nodes, no adjustment of hopping constants 
or wavefunctions amplitudes is needed. 

The fourth rule is that the amplitude at a given site 
of a localized state with zero energy can be renormalized 
without changing the energy of the state, if the hopping 
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(b) Rule II: add/remove nodes. 
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(c) Rule III: folding (d) Rule IV: renormalization 
(« = l/v/2). (E = 0). 



(e) Rule V: unfolding (a = 1/\/N). 


Figure 3. Set of origami rules for the construction of localized 
states in decorated lattices. The values of a can be obtained 
from simple tight-binding calculations, imposing the condi¬ 
tion that the state is still an eigenstate after applying the 
rule. 


constants to that site are renormalized as well (see Fig. 


3dl. 


The fifth rule describes the unfolding of a plaquette 
around a given site (see Fig. 3e). Multiple rotated copies 
of the original plaquette can be added around a site, pro¬ 
vided that the amplitude of the wavefunction on this site 
is adjusted, as well as the hopping constants around this 
site. 
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Figure 4. Localized states in the (a) Mielke lattice and (b) 
Tasaki lattice. All hopping constants are the same {tij = t), 
except those associated with curvy lines {tij = t/2). Note 
that the tight-binding Hamiltonian for the Mielke plaquette 
is symmetric in the exchange of sites A and B, but the full 
tight-binding Hamiltonian is not. This implies the Hamilto¬ 
nian eigenfunctions must have the same amplitude value (or 
opposite values) at sites A and B. In the case of the localized 
states, the value must be the same and the sites A and B are 
effectively one site. 


This set of rules justifies the existence of localized 
states in the Mielke and Tasaki lattices of Fig. In 
Fig. [I we exemplify the application of this set of rules 
starting from the localized state of the Lieb plaquette 
and ending at the localized state of the Tasaki lattice. 

We emphasize that these rules can be applied to con¬ 
struct localized states in systems of arbitrary dimension, 
from OD (a molecule) to 3D crystals, since the tight- 
binding bonds of Fig. may not be coplanar and the 
unfolding axis can have an arbitrary direction. Further¬ 
more, the spin degree of freedom of a tight-binding model 
may be interpreted as a sublattice index, and spin flip¬ 
ping terms can be interpreted as hopping terms between 
such sublattices. Localized states can also be created us¬ 
ing the above origami rules involving the two (up and 
down spin) sublattices and examples of such localized 
states are the edge states in ID models of topological in¬ 
sulators with spin flipping terms (see for example section 
3.3. of Ref. [321 • 

Another context where this set of rules applies is that 
of one-magnon localized states in frustrated quantum 
Heisenberg antiferromagnets. In these systems, these 
states generate peculiar behavior s uch a s magnetization 
plateaus around the saturation fleld ^^F^ I 

Let us now discuss how much of this method can be ap¬ 
plied in the U —>■ oo limit of the Hubbard model. We start 
by considering the Lieb lattice and then we generalize 
our conclusions to more complex decorated lattices. As 
discussed previously, the Lieb plaquette is a 8-site tight- 
binding ring. The eigenfunctions of the Hubbard ring in 
the strong coupling limit can be written as a tensorial 
product of the eigenfunctions of a tight-binding model of 
independent spinless fermions (holes) in the ring with L 
sites and the eigenfunctions of an Heisenberg model (with 


exchange constant J = /U) in a reduced chaiiJ^^Eni 

The ring of spinless fermions is threaded by a fictitious 
magnetic flux, (j) = Qs, generated by the spin configura¬ 
tions in the reduced Heisenberg chain (where Qs is the 
total spin momentum) and the eigenvalues to order t are 
given by 

E{ki,...,kN,^+Na) = -2t'^cos (ki - , (3) 

i=l 

where ki = {2 'k/L) ni, = 0,..., L — 1 are the momenta 
of the holes in the spinless ring. It is obvious that the 
rules for the construction of tight-binding localized states 
discussed above also apply to the case of one hole in a 
saturated ferromagnetic background. What one also con¬ 
cludes from the solution of the Hubbard ring in the strong 
coupling limit is that the same applies to the case of one 
hole moving in a Lieb plaquette with arbitrary spin con¬ 
figuration as long as the spin momentum is zero (since 
non-zero spin momentum destroys the time reversal sym¬ 
metry of the tight-binding model of spinless fermions). 

Does this apply to more complex plaquettes that share 
the rotation symmetry of the Lieb plaquette? Taking 
the example shown in Fig. one sees that in the non¬ 
interacting case, a localized state is present where the 
particle is confined to a ID path. This leads one to sug¬ 
gest that an equivalent localized state can be constructed 
for the hole moving in the spin background, if we impose 
a ( 7 s = 0 spin momentum for the spins configuration in 
the ID path. However one should note that in Fig. 
despite the electron probability density being finite only 
in the outer ring, when the electron is at the outer ring, 
it still hops to the center site, but summing over all the 
hopping possibilities from the sites at the outer ring to 
the center site, the result will be zero (destructive in¬ 
terference). In the case of the hole moving in the spin 
background, the hops of the hole from sites at the outer 
ring to the center site mix the spins at the outer ring 
and at the center. In order for one to have destructive 
interference at the center, the spin configuration must 
be a uniform linear combination of all possible permuta¬ 
tions of the set of spins (given the number of up spins 
and down spins)P^. The reason is the following: when an 
electron in the localized state of Fig. |^hops from a site 
of the outer ring to the center site, it interferes destruc¬ 
tively with the contributions of hoppings from the other 
sites of the outer ring. In the case of the hole, one has 
the different spin backgrounds and the hops of the hole 
from the outer ring to the center should not apparently 
interfere destructively. However, if one works with hole 
states such that, for a given number of up and down spins 
in the plaquette, the spin configuration is a uniform lin¬ 
ear combination of all possible permutations of the set of 
spins, then the spin configuration in the outer ring has 
Qs — 0 spin momentum, that is, one can freely perform 
circular shifts of the spin configuration in the outer ring. 
This implies that, assuming a spin configuration which 
is a uniform linear combination of all possible permuta- 
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Figure 5. Application of the rules presented in the text, starting from the localized state in the Lieb plaquette and ending at 
the localized state in the Tasaki lattice. Note that the energy of the localized state in the Tasaki lattice is determined by the 
central hoppings between sites with finite probability density. 


tions of the set of spins, a hole jump from a site A or B of 
the outer ring to the central site will generate the same 
final state, independently of the initial site being A or B. 
Therefore, we have the same localized state for one hole 
in the U ^ oo limit of the Hubbard model as for one 
electron when U = 0. This construction of a localized 
state is valid for any plaquette corresponding to a deco¬ 
rated lattice where a localized state of one tight-binding 
electron exists. 

Note that besides the localized state degeneracy asso¬ 
ciated with the choice of the lattice plaquette where the 
localized state sits, there is a huge degeneracy associated 
with the possible choices of number of up spins (or down 
spins) in the plaquette and in the rest of the lattice. This 
degeneracy is lifted by the Heisenberg corrections of order 

(U, as in the Hubbard rin^^^H^. 

Another important remark is that while the one- 
electron localized states of the Mielke and Tasaki lat¬ 
tices are the ground states of the respective tight-binding 
Hamiltonians, in the case of the U —>■ oo Hubbard model, 
the one-hole localized state is not the ground state, for 
the choice of relative hopping constants of Fig. How¬ 
ever, since the exact form of the hole probability density 
is known as well as the hole wavefunction phase (as in the 
case of the one-electron localized state), it is possible to 
tune the geometry, hopping constants and interactions 
in order to lower the energy of the hole localized state 
relatively to the other states, so that the energy of the 
localized state approaches the energy of the ground state. 

In conclusion, we have presented a simple set of rules 
for the construction of localized states of the Hubbard 
model in nearly arbitrary decorated geometries, in the 
tight-binding limit (U = 0), and in the strong-coupling 
limit ([/ —>■ oo). The first step in this method is the choice 
of a plaquette or a set of plaquettes with a higher symme¬ 
try than that of the whole lattice. In this plaquette, one 
has a localized state of the tight-binding Hamiltonian of 


the full lattice (this state has probability density nodes 
at the sites shared between the plaquette and the rest of 
the lattice). Using a simple set of rules, the tight-binding 
localized state in such plaquette can be divided, folded 
or unfolded to new plaquette geometries. We have shown 
that this set of rules can also be applied in the U —>■ oo 
limit of the Hubbard model, for the construction of local¬ 
ized states of one hole in the plaquette, assuming a spin 
configuration which is a uniform linear combination of all 
possible permutations of the set of spins in the plaque¬ 
tte. Note that in every other plaquette, one may place a 
localized hole, so localized hole states exist for hole dop¬ 
ing between zero and a value of the order of l/A’piaq (the 
value depends on the lattice geometry), where fVpiaq is 
the number of plaquette sites. 

This paper presents a unifying picture of construction 
of localized states, in tight-binding systems of arbitrary 
dimension (from OD to 3D), arbitrary geometry (includ¬ 
ing Mielke’s and Tasaki’s 2D geometries), without and 
with interactions (U = 0 or U = oo, extending in the 
latter case the filling intervals where localized states are 
known to occur). The existence of localized states due to 
spin flipping terms in tight-binding descriptions of topo¬ 
logical models, or the existence of one-magnon localized 
states in frustrated Heisenberg antiferromagnets, are two 
other contexts included in this unifying picture. 
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